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pilot's  reference  frame  and  adjoins  these  to  the  aircraft  attitude 
states  and  two  oscillating  wind  states. 
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Ad  st  rsc  t 

Although  the  accuracy  of  the  srufi  . riiij;  .v . 1 5 C 
aircraft  is  excellent  relative  to  other  weapon  systems, 
a detailed  error  analysis  reveals  that  the  pilot  is  by 
far  the  primary  source  of  remaining  error.  This  report 
concludes  that  the  accuracy  of  his  pointing  the  aircraft 
could  be  dramatically  improved  if  the  sight  line  were 
aimed  automatically.  The  sight  line  autopilot  (SLAl’J 
is  designed  to  meet  this  need.  The  autopilot  design 
accounts  for  the  linearized  motions  of  the  sight  line  as 
seen  in  the  pilot's  reference  frame  and  adjoins  these 
to  the  aircraft  attitude  states  and  two  oscillating 
wind  states. 

Using  optimal  regulator  theory,  the  aircraft  control 
gains  are  generated.  Using  these  gains,  extensive  computer 
simulations  were  run  to  check  the  autopilot's  control. 
Included  were  winds  and  sens  ir  noise.  The  results  show 
that  the  concept  of  a Sl.AP  will  signi  ficant  iy  improve 
the  capability  of  the  sidefiring  weapons  system  and  allow 
the  pilot  to  concentrate  on  aiming  the  weapons. 
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INTRODUCTION 


The  development  of  the  gunship  concept  for  air  warfare 
has  progressed  steadily  since  the  early  operation  of  the 
AC-47  "Dragonships"  through  the  initiation  of  the  AC- 130 
model  aircraft  into  the  program.  The  degree  of  sophistica- 
tion has  also  increased  from  an  area  coverage  device,  which 
required  little  accuracy,  to  a designated  target  device, 
requiring  a high  degree  of  accuracy.  This  accuracy  require- 
ment demanded  a further  investigation  of  error  sources  in 
the  weapon  system. 

The  U.S.  Air  Force  Academy,  as  part  of  the  program  of 
research  and  consulting,  has  been  steadily  supplying 
assistance  to  the  Gunship  Office  and  was  requested  to 
conduct  the  principal  investigation  of  error  sources  and 
their  possible  correction.  Early  research  determined 
that  the  pilot  was  a major  error  source.  This,  however, 
was  not  entirely  his  fault.  The  AC- 130  has  several 
slightly  damped  dutch-roll  characteristics  which  make  it 
very  difficult  to  hold  aircraft  attitudes  to  within 
fractions  of  a degree.  Vastly  complicating  the  picture 
is  the  fact  that  the  attack  attitude  is  a steady  turn  with 
the  guns  pointed  out  the  side  of  the  aircraft.  Thus, 
lateral  and  longitudinal  motions  are  not  separable. 
Furthermore,  the  required  attack  mode  is  to  fly  at  nearly 
constant  altitude  and  air  speed  on  any  given  firing  run. 
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The  requirements  to  maintain  the  necessary  flight 
condition  to  achieve  the  accuracy  are  so  severe  that  the  pilot 
would  apparently  benefit  from  a specialized  autopilot.  This 
autopilot  could  hold  the  weapons  on  target,  even  in  the  face 
of  winds,  using  a target  line.  The  important  thing  is  to  be 
able  to  fire  and  hit  the  target  --  not  to  fly  a circular 
pattern.  Thus,  with  this  autopilot,  the  target  should  always 
lie  along  the  target  line  which  is  only  slightly  affected  by 
fire  control  corrections  or  offsets  from  a constant  direction. 

A problem  overview  can  be  obtained  from  Figure  1.  It  can  be 
seen  that  this  is  a classic  regulator  problem.  The  principle 
complication  is  that  this  is  at  least  a twelve  state,  three 
input  system.  Because  of  the  turning  condition,  it  is  not  easy 
to  separate  the  modes.  Furthermore,  because  of  the  interest 
in  maintaining  a constant  target  line,  the  wind  is  a very 
significant  factor. 

Using  linearized  dynamics,  it  is  possible  to  cast  this 
problem  in  an  optimal  regulator  format  that  uses  a Riccati 
equation  to  establish  state  feedback  gains.  The  optimal  design 
has  been  verified  by  extensive  simulations  including  the  effect 
of  sensor  noise,  outside  unknown  constant  error  sources,  and 
limiters  on  the  control  surface  motion.  The  digital  program 
principally  used  in  this  investigation  is  the  Linear  Quadratic 
Loss  program  created  at  the  University  of  Florida,  by,  among 
others.  Captain  Charles  Fosha  (Department  of  Astronautics  and 
Computer  Science,  USAFA)  as  an  evolution  of  previous  ASP  and 
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GASP.  Extensive  modification  to  the  program,  notably  the 
technique  of  spectral  factorization  to  solve  the  matrix  Riccati 
equation,  was  accomplished  by  Captain  Fosha  in  support  of  the 
autopilot  design.  As  is  true  of  all  numerical  programs;  certain 
ruses,  described  later  in  this  paper,  were  employed  to 
avoid  numerical  difficulties. 
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<$E  * Elevation  Angle  between  Sight  Line  and  Target  Line 
5 * Azimuth  Angle  between  Sight  Line  and  Target  Line 

A 

X -Y  plane  is  EARTH 
0 0 

Figure  1 - Problem  Geometry 
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The  dynamic  model  equations  for  the  gunship  problem  can 
be  derived  in  two  parts.  These  two  parts  correspond  to 
(1)  the  motion  of  the  aircraft  with  respect  to  a given  air 
mass,  and  (2)  the  motion  of  a point  on  the  ground  with 
respect  to  an  observer  on  the  aircraft.  The  development  of 
the  linearized  equations  of  motion  for  an  aircraft  about  a 
nominal  flight  condition  are  well  known  and  shall  not  be 
repeated  here.  For  a complete  derivation  see  Dynamics  of 
Flight  by  Bernard  Etkin  (Ref). 

The  linearized  states  which  are  used  are: 
xi  = 6a  (perturbation  to  angle-of-attack- radians) 
x2  = 6q  (perturbation  to  pitch  rate- rad/sec) 
x 3 = 6B  (oerturbation  to  sideslip-radians) 
x4  * 6p  (perturbation  to  roll  rate-rad/sec) 
x5  = 6r  (perturbation  to  yaw  rate-rad/sec) 
x6  <5 (perturbation  to  roll  angle- radians) 
x7  60  (perturbation  to  pitch  angle-rad'  ans) 

The  perturbations  are  defined  as 
6a  = a - aQ 

where  a is  actual  angle-of-attack  and  aQ  is  the  nominal  angle 
of  attack.  Similar  definitions  hold  for  the  perturbations  to 
the  remaining  states.  Since  the  nominal  flight  condition  is 
a level  equilibrium  turn,  the  nominal  values  used  are  constant. 
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The  aircraft  equations  were  used  as  follows,  assuming 
that  nominals  sideslip  (go)  is  zero. 

mu  <=  (Pu0SCx  ♦ l/2puQSCx  _ mq0“0)u  * (l/2puQ2SCx  - mqQu0)a 


+(l/4puQScCx  - nwQuo)q  - (mg  cos  9 )0  + (mr0u0)g 

q 

•(pu.SC  + l/2pu  SC  ) u - l/2pu  2S(C  ♦ -££-  C )ag 
° * ° xu  g a 8uQb  q 

-(l/4puQScCx  )q’g 

q 

mcx0u  + (muQ  - l/4p  uQScCz . )a  = (puQSCz  + l/2pUQSCz  ♦ mqQ)u 

+ (l/2puo2SCZa)a  ♦ (l/4pu0ScCz  + mu0)q  - (mg  cos  4>Q  si.i  9q)0 

-(•nu0P0H  - (mg  sin  4>Q  cos  60)<t>  + (l/2puo2SCz  )6e 

6e 

-(PuoSCz  + l/2puoSCz  )u  - l/2pu02S(C  ♦ C )ag 


a 8 b 4q 
o 


l/4puoScCzq)q*g 


- (1/ 4pUQSc  Cm  )d  + I 
d 


U m . yy  • o m 

r 2 " ‘2 

-j(l/2puo  ScCm  )a  ♦ (l/4pu  Sc  C 

••  /M  W 


yyq  " (^aoScCm  + l/Zp^ScC^u 


l)q  + frzz  ’ lxx)  To  ' 2IxzP^P: 


(Izz  ' + 2Ixzro]r  ♦(1/2pu02scCp  )5e 


-(pu0ScCm  * l/2p»0ScC^)ug  - l/2puo2Sc(Cma  ♦ -K_  C^ag 
-Cl/4puoSc  2CmQ)q«g 


6 - (cos  4» 0 ) q - (sin  <$>0)r  - (qQ  sin  4>o  + rQ  cos  4»q ) 4. 


i 


+ 


IXXP  ' XxzT  “ kSuobC£)u  + t^yy  ' 1 zz^ ro  + Txzpo]  q 

+ (l/4pu0Sb2Cjlp  + Ix2q0-)p  + [1/4pUoSb  CJt  r * ^Izz*Iyy^qo]  r 

2 

+ (l/2puo2SbC£  )3  ♦ (l/2puo^bC£  )6r  + (l/2puo  SbC^  )«a 
^ r ^ 

-(pSu0bCt)ug  - (l/<pu0Sb2Clp)pg  - (l/%,u02Sb2C1I.)r'B 

-l/2pu02Sb(C  - j-  Clr)6g 

oil 

o 

-IX2P  + IZ2r  = (pSuobCn)u  + (cixx  ■ Iyy5po  ' IxzTcl1  q 

♦ [l/4puoSb2Cnp  - (Iyy  - Ixx)  qj  P + d/4P%Sb2Cnr  - Ixzqo)r 

-Cl/2puo2SbCnp)3  + (l/2pu02SbCn  )6r  ♦ (l/2puo2SbC^  )« a 

4"  a 

-(pSu0bCn)ug  - (l/4pu0Sb2Cn^)Pg  -(l/4pu02Sb2Cnr)r’g 
-l/2pu  2Sb(C ne--5—  C )6g 


muoe  = (Pu0SCy  - mro  + mpooto)u  +(mpouo)<x  * (mg  sin  Sin  6o)9 

+ (l/4pu  SbC  ♦ ma0a0)p  + (l/4pu0SbCy  - mu^r 
XP  r 

2 

+ Cl/?  ,u02SCy  )B  + (nig  COS  cos  eQ)4>  +(l/2puQ  SCy^  - (pu0SCy)ug 
3 r 

- (l/4pu0SbCy^)Pg  - (l/4pu0SbCyr)r'g  - l/2pu02S(Cy^  - -1-Cy^)3g 

i = (sin  $Qtan  6Q)q  + [(qQsin  4>0  + rQ  cos  <t>0)sec2  6 J 6 + p 

n (cos  4>0  tan  6Q)r  ♦ [(qQcos  <f>0  - r0  sin  4>0)  tan  6Q]  * 


The  aircraft  states  which  are  not  used  are  azimuth  (4>)» 

(which  is  cyclic)  and  velocity  variation  (6u)  which  is  main- 
tained at  zero  by  the  copilot  or  a separate  control  device.* 

POSITION  DYNAMIC  EQUATIONS 

The  gunship,  during  its  attack  maneuver,  attempts  to  main- 
tain the  target  in  a very  nearly  fixed  position  in  the  air- 
craft coordinate  frame.  This  position  is  on  the  left  beam 
and  below  the  wing  line.  Thus,  the  aircraft  flies  a modified 
pylon  turn  about  the  target.  For  this  reason  it  is  natural 
to  select  an  aircraft  oriented  coordinate  frame  in  order  to 
define  the  translational  position  of  the  aircraft.  In  other 
words,  the  autopilot  will  attempt  to  position  the  target  out 
the  left  wing  and  down,  using  the  coordinates  Range  (R) , 

Elevation  (E) , and  Azimuth  (A)  which  are  defined  in  Figure  2.** 

A second  independent  set  of  translational  coordinates  is  alti- 
tude (H) , with  E and  A as  previously  defined.  The  change  of 
reference  frames  from  range  to  altitude  is  a simple  transfor- 
mation and  is  accomplished  later  in  this  section.  The  three 
additional  states  are  designated  as  follows: 

x8  * R (range- feet) 

x*  » E (elevation  angle  of  sight  line-radians) 

xio=  A (azimuth  angle  of  sight  1 ine - radians) 

A forcing  function  due  to  wind  will  be  included  in  the  derivation. 

*This  state  will  probably  be  added  to  the  problem  along  with 
a throttle  control  variable  at  a later  time. 

**Note  the  nominal  elevation  angle  is  negative,  nominal  azimuth 
is  zero  and  nominal  range  is  the  slant  range  to  the  target. 
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The  time  derivative  of  the  radius  vector  in  the 
earth  frame  consists  of  the  velocity  of  the  aircraft  in 
the  air  mass  and  the  motion  of  the  air  mass.  Expressed 
in  the  aircraft  oriented  frame,  the  velocity  of  the  air- 
craft with  respect  to  a point  in  the  air  mass  is 


u cos  a cos  8 
u cos  a sin  8 
u sin  a cos  8 


where  u is  the  total  velocity  of  the  aircraft  in  feet  per 
second. 

If  the  target  is  considered  a stationary  point  on  the 
ground  and  the  winds  are  constant,  the  velocity  of  the 
air  mass  with  respect  to  the  target  is 


-Vw  cos  ip 
Vw  sin 


where  Vw  is  the  horizontal  wind  component,  is  the  vertical 
wind  component  and  ip  is  the  angle  in  the  horizontal  plane 
from  the  wind  direction  to  the  projection  of  the  aircraft 
x-direction  (nose).*  So,  the  velocity  of  the  aircraft 
with  respect  to  the  target  and  the  velocity  of  the  air- 


*<P  = aircraft  heading  - wind  direction  s 


a 
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craft  with  respect  to  the  air  mass  can  be  expressed  as: 


u 

cos 

a 

cose 

■ 1 

0 

0 

- 

cose 

0 

- s in0 

Tvw 

cos  4; 

u 

cos 

a 

sine 

+ 

0 

cos 

sin 

♦ 

0 

1 

0 

’vw 

sin  t 

_u 

sin 

a 

cose. 

0 

* 

-sin$ 

cos 

♦ 

sin0 

0 

COS0 

> 

> 

<— 
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The  0 and  $ transformations  are  required  to  take  the 
wind  velocity  from  the  horizontal  plane  to  the  aircraft 
(x,y)  plane  as  illustrated  in  Figure  3.  Since  this  represents 
the  velocity  of  the  aircraft  with  respect  to  the  target 
as  measured  in  the  earth's  frame,  the  velocity  of  the  target 
with  respect  to  the  aircraft  requires  only  a sign  change. 

The  result  is 


! 

r-u 

COS 

a 

cos 

en 

“i 

0 

0 

E* 

R = 

-u 

COS 

a 

sin 

e 

+ 

0 

cos  <J> 

sin 

<P 

J. 

-u 

L 

sin 

a 

cos 

0 

-sin$ 

cos 

♦_ 

“ cos 

G 

0 

-sin 

e' 

"vw 

cos 

* ~ 

0 

1 

0 

*vw 

sin 

* 

sin 

0 

0 

cos 

e 

1 

< 

< 

R 

where  K is  the  radius  vector  from  the  aircraft  to  the 

target  as  shown  in  Figure  2. 

The  "Coriolis  Law"  is  now  used  to  find  the  relation- 
E R 

ship  between  R and  the  time  derivative  of  "K  in  the  aircraft 
K 

* 

( ) 

X - Vector  X;  expressed  in  [ ] corrdinates,  with  time 
[ ] derivative  in  ()  coordinates. 
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frame 


A 

R 

X 


E A R 

R * R + u)  x X 
X X A/E 


(2) 


w^E  is  the  angular  rate  of  the  aircraft  with  respect 
to  the  earth. 

R 

Referring  to  Figure  2,  X can  be  represented  as 


R 

X 


R cos  E sin  A 
-R  cos  E cos  A 
-R  sin  E 


(3) 


The  aircraft  angular  rate  vector 


-A/E 


P 

q 

r 


can  be  expressed  as  the  matrix 


0 

r 

q 


~r 

0 

P 


q 

-p 

0 


(wA/Ex) 


(4) 
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to  evaluate  the  cross  product  terra  in  matrix  muliplication 
form;  i .e.  , 

D 

^A/E  x X *s  reP^ace^  with 


‘ R 

cos 

E 

sin 

-R 

cos 

E 

cos 

-R 

u 

sin 

E 

Now,  to  fill  in  for  R,  the  time  derivative  of  Equation 

/ » i * __  x l . ^ : r*.  c M 


(3)  as  run  in  the 


aircraft  frame  is  taken 


• * • 

R cos  E sin  A - R sin  E sin  AE  + R cos  E cos  AA 

* • 

= -R  cos  E cos  A ♦ R sin  E cos  AE  ♦ R cos  E sin  AA 
• • 

-R  sin  E - R cos  EE 

Equation  (2)  is  now  represented  as 


■u  cos  a cos  Bl  1 0 0 

•u  cos  a sin  B ♦ 0 cos$  sin<|> 


u sin  a cos  81  |J)  -sin$cos$j  ^sin® 


cos6  0 -sin6  Vw  cos 
01  0 -Vw  sin  <l> 


R cos  E sin  A - R sin  E sin  AE  + R cos  E cos  AA 
• • • 
-R  cos  E cos  A + R sin  E cos  AE  ♦ R cos  E sin  AA 
« • 

-R  sin  E - R cos  EE 


cosBJ  |_V 
(6) 


R cos  E sin  A 
-R  cos  E cos  A 
-R  sin  E 


Equation  (6)  can  be  solved  for  the  line-of-sight 


states 


and  linearized  about  given  nominal  conditions.  Assumptions 
regarding  the  nominal  condition  are: 

1)  small  angle  approximations  used  for  a e , 

<5E,  <5A  oo 

2)  B0  - 0,  Ao  = 0 


The  result  is  shown  in  Equation  (8) 
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THE  WIND  DYNAMICAL  EQUATIONS 

From  the  seven  air  mass  states  and  the  three  line- 
of-sight  states,  we  define  two  states  which  consider  the 
wind  effects : 

xi  1 * Vw  cos  \p  (wind  along  the  A/C  heading) 

Xj 2 * sin  1 p (cross  wind) 

Taking  the  derivatives  of  the  wind  states  and  linearizing 

xn  “ Vw  sin  *xi  2^  = -xi2^0 

X1  2 ■ Vw  COS  tplp  = X1  llp^xi  li0 

CONVERSION  TO  H.  A.  E 

There  is  a second  set  of  position  states  which  is 
possibly  more  useful,  depending  upon  the  availability  of 
sensors,  since  we  are  ultimately  concerned  with  state  feed- 
back. This  second  set  uses  altitude  H,  and  the  two  angles 
E and  A.  To  derive  the  second  set,  two  relations,  one 
static  and  one  dynamic,  are  established  for  altitude. 

The  static  relation  is  obtained  by  transforming  the 
vector  through  the  angles  0 and  <J>  to  the  earth  frame 
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Altitude  is  found  in  the  third  component: 

-R  cos  E sin  A sin  6 - R cos  E cos  A sin  ♦ cos  6- 
R sin  E cos  $ cos  6 

Using  small  angle  approximations  for  6 and  A and  a 
trigonometric  identity  for  sin  (<|>+E),  we  obtain 

H * -R  sinO+E)  - R8A  cos  E 

The  dynamic  relationship  is  similarly  obtained  by 
transforming  the  aircraft  velocity  in  the  air  mass 
through  the  angles  $ and  0 and  determining  the  third 
component.  The  result  is: 

• 

H * u cos  a cos  6 sin  0-u  cos  a sin  g sin  $ cos  0- 
u sin  a cos  g cos  cos  0 

Using  the  small  angle  approximations  for  a,  3 and  <J> , 
results  in 

H * u9  - ug  sin  <J>  - ua  cos  d> 

Equation  (8)  can  now  be  replaced  with  Equation  (9),  which 
is  on  the  following  page. 

The  seven  aircraft  states,  the  three  position  states, 
and  the  two  wind  states  form  the  dynamic  plant  for  the 
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problem.  The  control  elements  chosen  are  the  aircraft 
control  surfaces.  That  is: 

U]  = 6Ru;  Rudder 

U2  = 6Ae ; Aileron 

U3  = SE^;  Elevator 

Conspicuously  left  out  here  is  a throttle  control; 
however,  recall  that  for  other  reasons,  the  throttle  or 
(6u)  is  maintained  at  zero. 

The  system  we  have  defined  here  is  of  twelve  coupled 
states  with  three  independent  inputs.  This  forms  a 
linear  dynamic  plant.  The  next  section  details  the 
problem  formulation  in  the  classic  linear  optimization 
method,  and  explains  some  details  of  the  optimization 
method . 


r*,  r . X-  1 nni  ; 
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FEEDBACK  CONTROL  SYNTHESIS 

The  problem  being  addressed  is  that  of  a linear 
regulator.  The  solution  for  quadratic  penalties,  using 
optimal  control  theory  is  explained  in  reference  2, 
Bryson  and  Ho.  The  results  will  be  briefly  reviewed 
here  to  establish  notation.  Given  the  system 


with 


x = Fx  + Gu 

x and  nxl  state  vector 
u an  mxl  control  vector 
F an  nxn  system  matrix 
G an  nxm  control  distribution  matrix 
minimize  the  Performance  Index  (J) 
t, 


J = hC  (xT  Ax  + uTBu 

^t 


Bu)dt 


where  A is  an  nxn  positive-semidefinite  matrix 
and  B is  an  mxm  positive  definite  matrix 
The  steady  state  solution;  i.e.,  for  large  t£,  is 
found  by  solving  the  Riccati  equation  "backwards"  in  time: 

S = -SF  - FTS  + SGB" 1 GTS  - A (10) 

until  it  reaches  a stabilized  value  of  S(a  constant) . 

This  S is  then  used  to  specify  the  feedback  matrix  C: 
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u * Cx 

c * -b‘‘gts 

So  the  optimum  system  "closed  loop"  equation  is: 
x = rF  - GB_1GTS]x 

This  theoretical  problem  is  solved  at  USAFA  using 
the  Linear  Quadratic  Loss  program.  The  solution  to  the 
steady  state  version  of  (8)  is  performed  algebraically 
using  Potter's  method.  This  depends  on  solving  for  the 
eigenvalues  and  eigenvectors  of  a 2n  x 2n  matrix.  Un- 
fortunately, repeated  roots  give  rise  to  numerical  diffi- 
culties, as  do  undamped  roots.  The  follouing  problems 
were  encountered  as  a result. 

a.  Regulating  to  non- zero  values.  Each  state  which 
is  called  on  to  regulate  to  a non- zero  value  can  be  handled 
by  adding  an  additional  integrator  which  has  an  initial 
condition  equal  to  the  desired  offset  value,  and  no  input. 
The  state  weighting  matrix  (A)  is  then  expanded  to  include 
the  additional  state,  and  the  quadratic  form  (A)  is  adjusted 
tc  weight  the  difference  between  actual  and  desired 
squared . 

This  led  to  two  numerical  problems.  There  are  three 
off-nominal  commands  that  are  useful:  those  for  azimuth 

(A-state  10),  for  altitude  (H-state  8)  and  for  elevation 
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(E-state  11).  The  three  integrators  represent  three  I 

> 

repeated  zero  roots  and  prevented  a numerical  solution.  ■ 

i 

This  was  solved,  however,  by  lightly  damping  each  of  the  ! 

i 

- I 

integrators,  (i.e.,  time  constants  of  10*  seconds)  | 

i 

slightly  differently.  There  was  no  appreciable  change 
in  state  during  the  duration  of  the  problem.  However,  ! 

the  eigenvectors  could  now  be  found  using  a double- 
precision routine  on  the  computer.  The  second  difficulty 
was  that  numerical  rounding  would  sometimes  cause  the 
weighing  matrix  A to  appear  indefinite  to  the  computer. 

This  was  solved  by  allowing  the  diagonal  terms  to  slightly 
dominate  the  off-diagonal  ones. 

b.  The  Wind  Oscillator.  Because  the  eigenvalues 

of  the  wind  oscillator  (states  11  and  12)  are  purely  imagi- 
nary, the  Potter  method  creates  a repeated  pair.  This 
is  because  all  system  eigenvalues  and  a corresponding  set 
which  is  reflected  about  the  imaginary  axis  are  present. 

This  problem  is  solved  by  adding  a very  slight  amount  of  ! 

damping  to  the  wind  oscillator. 

c.  Hydraulic  Servo  Lags.  To  validate  the  performance 

i 

of  the  autopilot  design,  servo  lags  estimated  to  be  .1 

second  were  added  to  each  actuator.  Digital  difficulties 

due  to  the  repeated  roots  were  quickly  resolved  by  adding 

a small  separation  between  the  three  time  constants.  , 

i 

\ 

> 

r 
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The  selection  of  the  matrices  A and  B was  initially 
made  on  the  basis  of  the  Bryson  rule  (reference  2,  page 
149).  That  is,  weights  on  A,  E and  H which  were  inversely 
proportional  to  the  acceptable  maxima;  and  weights  on 
control  which  are  inversely  proportional  to  the  maximum 
desired  control  deflections.  The  resulting  system  tended 
to  respond  so  quickly  that  rigid  body  assumptions  were 
suspect.  To  solve  this,  weights  were  added  to  the  body 
angular  rates  until  the  digital  simulation  demonstrated 
reasonable  behavior  for  the  assumed  range  of  initial 
conditions . 

RESULTS 

In  order  to  verify  theoretically  the  control  solution 
obtained  by  solving  the  linear  regulator  problem,  three 
cases  were  simulated.  The  cases  presented  are  in  order 
of  increasing  difficulty.  The  cases  were  simulated  for 
100  seconds,  or  approximately  two  orbits. 

Case  1 - Error  in  Position  Given  No  Wind  and  No  Exterior 
Disturbance 

The  first  case  is  presented  as  a basis  for  comparison 
against  the  more  difficult  second  and  third  cases.  The 
initial  conditions  are  as  follows: 

6h  = 100  ft;  6A  - 5E  = 3° 
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The  results  in  this  case  are  given  in  Figure  4 (a  5 b). 

The  influence  of  weighting  the  p,  q and  r attitude  states 
is  immediately  apparent  as  the  position  errors  slowly 
damp  towards  zero.  Without  this  weighting  the  practicality 
of  the  controller  would  be  negligible  due  to  real  control 
saturation  and  excessive  "g"  forces  impacting  the  crew  of 
the  aircraft.  The  control  surface  deviation  is  given  in 
the  lower  graph  anl  is  well  within  practical  aircraft 
1 imits . 

Case  2 - Error  in  Position  Given  a Constant  Wind  But  No 
~ External  Disturbance 

Although  not  shown,  it  is  logical  that  the  path  flown 
wilr  not  be  wholly  circular,  but  rather  an  ellipsoid.  Case 
2 illustrates  the  capability  of  the  autopilot  to  maintain 
a near  zero  error  with  wind.  The  initial  conditions  are  as 
follows : 

wind  = 40  ft/sec  head  on 

6h  - 100  feet;  <5A  = <5E  = 3°  „ 


Presently,  in  the  constant  wind  condition,  the  pilot 
cannot  maintain  contact  with  the  target  for  greater  than 
one-half  a turn.  The  results  of  the  theoretical  simulation 
clearly  show  (Figure  5 (a  5 b))  that  contact  is  possible 
throughout  the  two-orbit  period.  The  errors  are  small  in 
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magnitude,  though  cyclic,  and  will  have  only  a minor  effect 
on  the  accuracy  of  the  weapons  system.  The  largest  relative 
error  occurs  in  altitude.  This,  however,  has  the  least 
effect  on  the  weapon  delivery. 

The  control  surface  deviations  are  still  small, 
however,  considerably  larger  than  those  in  Case  1.  It  is 
notable,  however,  that  there  are  no  rapid  changes  in  control 
surfaces,  even  at  the  start  and  that  the  largest  deviation 
remains  well  within  the  aircraft  operating  range. 

Case  3 - Error  in  Position  Given  a Gusting  Wind  and  External 
Disturbance  on  All  System  States 

The  noise  introduced  in  this  case  is  entered  as 
unknown  error  sources,  assuming  that  the  measurements 
remain  known.  The  nominal  errors  used  for  initial  conditions 
are  as  follows: 

wind  = 40  ft/sec  head  on 
<5h  = 100  ft;  6A  = SE  = 3° 

This  case  is  illustrated  by  Figures  6 (a  5 b).  The  path 
flown  by  the  aircraft  is  again  an  ellipsoid,  with  the  added 
difficulty  of  the  external  unknown  errors.  Once  again,  the 
control  gains  are  capable  of  maintaining  target  contact 
throughout  the  two-orbit  period.  The  accuracy,  as  expected, 
is  impaired  by  the  external  errors;  however,  the  number 
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The  sight  line  autopilot  is  part  of  a funded  Air  Force 
development  program.  A contract  has  been  let  to  Minneapolis 
Honeywell  to  determine  the  exact  details  based  on  the 
studies  performed  at  USAFA.  It  is  clear  that  all  state 
variables  should  be  instrumented.  To  accommodate  different 
gain  values  for  different  flying  configurations  a digital 
computer  is  very  desirable,  and  probably  will  be  part  of 
the  test  aircraft  configurations.  The  SLAP  will  probably 
parallel  existing  flight  hardware  and  consequently  offer  a 
redundant  control  method  in  the  event  of  battle  damage. 
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CONCLUSIONS 

The  complexity  of  the  autopilot  design,  and  the 
linearization  of  the  dynamic  equations  about  a nominal 
flight  condition  made  the  problem  ideal  for  optimal  control 
theory.  With  the  necessity  for  observing  all  states,  a 
Kalman  filter  becomes  a possibility. 

The  results  clearly  show  that  theoretically  the  addition 
of  the  sight  line  autopilot  will  greatly  enhance  the  accuracy 
of  the  gunship  operation.  The  results  do  not  show  the 
increase  in  pilot  potential  due  to  relieving  him  of  the 
targeting  problem.  The  resulting  controller  will  be  able  to 
theoretically  control  the  side-firing  aircraft  sight  line  to 
an  accuracy  of  better  than  one  mil.  This  development  should 
significantly  advance  the  capability  of  weapons  systems  of 
this  type. 
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